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Abstract 

In the framework of renormalization- group improved cosmologies, we use the Noether symmetry 
approach to get exact and general integration of the matter-dominated cosmological equations. 
This is performed by using an expression of A = A(G) determined by the method itself. We 
also work out a comparison between such a model and the concordance ACDM model as to the 
magnitude-redshift relationship, hence showing that no appreciable differences occur. 
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I. INTRODUCTION 



The study of cosmological dynamics has been recently performed by analyzing strong 
"renormalization group (RG) induced" quantum effects. These are supposed to drive the 
(dimensionless) cosmological "constant" X(k) and Newton "constant" g(k) from an ultra- 
violet attractive ffxedpoint. This approach acts within the framework of the effective 
average action 1, [2], 3|, finding a fixed point in the Einstein-Hilbert truncation of the- 
ory space 0, [5I, 6] and in the higher-derivative generalization Q]. The existence of such a 
non-Gaussian ultraviolet fixed point in the exact theory implies its nonperturbative renor- 
malizability 0, ?j, Isl, El ID, U, 12]. This RG-improved framework describes gravity at a 



typical distance scale £ = k^ 1 , introducing an effective average action r^fgr^] for Euclidean 
quantum gravity [4j , and deriving an exact functional RG equation for the fc-dependence of 
rV This context is usually referred to as quantum Einstein gravity, and makes it possible 
to find an explicit fc-dependence of the running Newton term G(k) and the running cosmo- 
logical term A(k). This is interesting for an understanding of the Planck era immediately 

nnn 

after the big bang as well as the structure of black hole singularity |13l . 114 , 1151 ] . 

In order to obtain the RG-improved Einstein equations for a homogeneous and isotropic 
universe, one can identify k with the inverse of cosmological time, k oc 1/t Q, Q], and 
a dynamical evolution for G and A(k) induced by their RG running can thus be derived. 
An Arnowitt-Deser— Misner (ADM) formulation has also been presented jl^]; it builds a 
modified action functional which reduces to the Einstein-Hilbert action when G is constant. 
A power-law growth of the scale factor can then be obtained for pure gravity and for a 
massless <£> 4 theory in a homogeneous and isotropic space-time, in agreement with what is 
already known on fixed-point cosmology. On the other hand, still within the framework of 
the ADM formulation, in 18j] we have proposed solutions for the pure gravity case derived 



by means of the so-called Noether Symmetry Approach. This is a method which implements 
a variable transformation that usually makes it possible to find exact and general solutions 
of the cosmological equations [13, |20[]. The solutions found in Ref. [18] predict that the 
Universe is in an accelerated stage, hence mimicking inflation without introducing a scalar 
field in the cosmic content. 

There are also some indications 21] that quantum Einstein gravity, because of its in- 



herent infrared divergences, is subject to strong renormalization effects even at very large 
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distances. In cosmology, such effects might be relevant for the Universe at late times, since 
they lead to a dynamical relaxation of A and, probably, may solve the cosmological constant 
problem 21|. As a matter of fact, the late accelerated expansion of the Universe can be 



viewed as a renormalization group evolution near a non-Gaussian infrared fixed point [16], 
where G and A become running quantities at some late time. The sharp transition between 
standard FLRW cosmology and accelerated RG driven expansion is supposed to occur at 
some time (the fixed point being reached exactly at the transition), but is a strong simpli- 
fying assumption; however, some agreement can be found between this kind of model and 
Snla observations 22fl. In what follows, we would like to present some exact solutions of the 

n 

improved Hamiltonian cosmology described in [17J], without any assumption on fixed-points 
structure, discussing the possible physical viability of the resulting model. Interestingly, 
it is possible to think that the kind of 'post-inflationary' accelerated stage we discover in 



our considerations could explain what is observed by means of Snla data |23l . |24|, |25|, |26 |. 
We stress again that we do not assume a priori what functional relation exists between A 
and G, but we merely postulate its occurrence. Some solutions compatible with having 
AG = constant were found in the ultraviolet regime in 171 ]. whereas we here find a different 
renormalizat ion-group trajectory. 

Here, we investigate the matter-dominated universe, limiting our analysis to cosmology 
in the dust case. Our aim is to derive the main parameters that make it possible to compare 
our model with the concordance ACDM model as to the magnitude-redshift relationship, 
hence showing that no appreciable differences indeed result when one wants to examine 
Type-la supernovae data. In section 2 we describe the Lagrangian formulation for the 
Lagrangian adopted to derive the RG-improved Einstein cosmological equations with the 
ordinary matter energy-momentum tensor and then find the Noether symmetry. Section 3 
studies how this gives rise to exact and general solutions, while in section 4 we discuss the 
situation with K, = 0, comparing our model and the ACDM model. In section 5 we draw 
conclusions. 



II. NOETHER SYMMETRY 



First of all, let us consider the approach outlined in Ref. 



171 ] and there applied to models 



of gravity with variable G and A. From the modified action functional, an Einstein-Hilbert 



3 



action can in fact be obtained when G and A are such. As said, in a homogeneous and 
isotropic universe, this leads to power-law behaviours of the cosmic scale factor a = a(t) 
for both pure gravity and a massless y? 4 theory, in agreement with results from fixed-point 
cosmology, once one adopts the constraint GA = const. On the other hand, it is known that 
an independent dynamical G is equivalent to metric-scalar gravity already at classical level 



271 . |28| . Taking both G and A as varying independently with position and time can drive to 



pathological situations, and it is possible to show that we must therefore assume A = A(G) 
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Here, we want to investigate the matter-dominated case in homogeneous and isotropic 
cosmology (with a si gna ture — , +, +, + for the metric, lapse function N = 1 and shift vector 
N l = 0). As in Ref. 17] (but see also Ref. 0), let us start from the Lagrangian 



8ttG 



1 G 2 

—3ad 2 + 3/Ca — a 3 A + -^a- 3 -^ 

2 G 



(2.1) 



where G = G(t), A = A(G(t)), JC = —1,0, 1 (for open, spatially flat and closed universes, 
respectively), while dots indicate time derivatives and /i is a nonzero interaction parameter 



introduced in Ref. [17( and also used in Ref. [18] . We have now inserted also the matter 
contribution by means of L m = -Da" 3 ' 7_1 \ with 1 < 7 < 2 (being, in the most relevant 
cases, 7 = 1 for dust and 7 = 4/3 for radiation) and D a suitable integration constant 
connected to the matter content. From Eq. (12. ip we find the second-order Euler-Lagrange 
equations for a and G 

" 2 /C A aG fiG 2 



a a 
a + 2a^ 



3 G 



2a 2 



aG AG 2 



4vrG(7 - l)Da-^ = 



G 



2 r G 



M 2 



QIC 



,dA 



/iG - -/i— + 3fi-G + - -6— + — - 2A + 2G— 



We have also to consider the Hamiltonian constraint 171 ] 



dG 



K, A [iG 



87TG 



^ + ^-4-^~^Da-^ = 



6G 2 



(2.2) 
(2.3) 

(2.4) 



which can indeed be seen as equivalent to the following constraint on the energy function 
associated with L 0,0,0: 



E, 



dL. dL- 
—a H -G 

da dG 



L = 0. 



(2.5) 
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In what follows we want to discuss the dust case, so choosing 7 = 1; this involves a 
zero pressure p m and an energy density p m = Da~ 3 . The Lagrangian L takes therefore the 
simplified form 



8ttG 



I G 2 
Sad 2 + 3/Ca — a 3 A + -/ia 3 



— D 



r ~ G2 , -, (2.6) 

so that the matter term is now just a constant and has no effect on the equations of motion 
with respect to the pure gravity case, but it is nonetheless important, since it always occurs 
in the constraint equation ( 12.51) . We have already solved the system of equations of motion 
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in pure gravity [18|, using the Noether Symmetry Approach [19J, [20J. Therefore, let us again 
consider the Lagrangian L as a p oint-like Lag rang ian, function of the variables a and G, 



and their first derivatives 



18 
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20J. In Ref. 18| . we have shown that, with a consistent 



choice of the function A = A(G), a Noether symmetry indeed exists for the pure gravity 
Lagrangian. A key point in our considerations is now given by the observation that, in the 
matter-dominated case, the procedure to deduce the symmetry is exactly the same as in the 
pure gravity situation, the only substantial difference being in Eq. (12.51) (which is in fact 
equivalent, now, to Eq. (12.41) ). therefore acting as a slightly different (but very important) 
constraint on the integration constants involved by the solution method. This means that 



we find exactly the same Noether symmetry as in 18 



so that we can here use the same 



transformations introduced there, and write the expressions for a = a(t) and G = G(t) 
obtained there, but taking care to consider, now, the updated energy constraint. 
More in detail, let us consider the vector field 



(2.7) 



with a = a(a,G) and (3 = (3(a,G) generic C l functions, and a = da/dt = (da/da)d + 
(da/dG)G, $ = d/3/dt = (d/3/da)a + (d(3/dG)G. As in Ref. |ia|, the condition 



CxL = 



(2.8) 



{CxL being the Lie derivative of L along X) corresponds to a set of equations for a 



a(a, G), (3 — (3{a, G) and A = A(G), which can be solved by 
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18| 



a(a, G) = as-tJG 



j-i 



P(a,Gt) 



3-2J 

The J parameter is an arbitrary constant 7^ 1,3/2 and such that 

M=|(3-2J)V0,|. 



a 3-2J G 



(2.9) 



(2.10) 



(In what follows we prefer to use J instead of /i so as to simplify the resulting expressions.) 
For consistency, we still get both 

J/C = (2.11) 



and 



2(W)A + G ^=0, 



(2.12) 



which splits our considerations into two separate branches. 
Thus, neglecting unnecessary integration constants, we have 



X, 



X|vj^o = a^G J ~ l — + aG J - 2 



.j-ijL 

da 



J . , , , N 3(J-i) . 

a 3 ~ 2J Ga + (J - l)a 3 - 2J G 



3-2J 



d_ 

da 



3 3(j-i) 7 d 

a a- 2 - 7 G — 1 

3 - 2 J dG 3 - 2 J 



a^-^G 



j-i 



3(J 



3-2J 



■Gd + JaG 



A = A(G) = WG 2(J_1) , 



A, (2.13) 
<9G 

(2.14) 



for J ^ 0, 1, 3/2 (=>■ any /x 7^ 0, 2/3) and K, = 0. It turns out that Xj — > Xq (as well as for 
the expression of A = A(G)) for J — > 0, but the situation with J = (=>• \i = 6) has to be 
nonetheless treated separately, since it is not a simple subcase except when /C = 0. 



III. SOLUTIONS FROM NEW COORDINATES AND LAGRANGIAN 

There exists a change of coordinates {a, G} — > {w,t> }, such that one of them (say w, for 
example) is cyclic for the Lagrangian L, and the transformed Lagrangian produces exactly 
and generally solvable equations. Solving the system of equations ixdu = 1 and ixdv = 
(where ixdu and ixdv are the contractions between the vector field X and the differential 



forms du and dv, respectively 
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20j), we get 

u = u(a,G) =na^G m , v = v(a, G) = In (aG~ 2nm ) . (3.1) 
Here, we have defined 

3 — 2 J 1 

nEn(J)E 6(M' m " m(J) " 1 " J= 2(3^1)' (3 ' 2) 

which are well defined and non- vanishing, being J 7^ 1,3/2. Eq. (I2.14p thus becomes 

A = A(G) = WG^ . (3.3) 
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In the case J = (i.e., n = 1/2 and m = 1) and any )C, we then get 



u = u(a, G) = —aG , 

v = v(a,G) = In (aG -1 ) . 

The inversion of Eq. (13. ip gives 

a = a(u,v) = n~ n exp ^— J-u™ , 

G = G(u,v) -- 



1 / v 
— exp I it 



l/(2m) 



which implies that 



a = a(u,v) = \/2uexp ^ 



when J = 0. Here, we set 



G = G?(u,u) 



2m exp 



-) 

2J' 

v 



Qn 2 



(1 - 3n) 2 ' 

The substitution of these functions a = a(u, v) and G = G(u, v) into L yields 

(m + 1 \ 
— f 1 — exp (3w) — 87rD 

for J ^ (and /C = 0), which becomes 

L' Q = —6 exp {2v)uv + 3/C exp (v) — W exp (3u) — 8tiD 



(3.4) 
(3.5) 

(3.6) 
(3.7) 

(3.8) 
(3.9) 

(3.10) 
(3.11) 

(3.12) 



for J = and any K,. By construction, in both cases u is cyclic for V , and there exists a 
non- vanishing constant of motion S = —dL'/dii associated to L', which helps to solve the 
equations deduced from it. We in fact find 

m+1 



6 exp 



2nm 



v I v , Tiq = 6 exp (2v)v , 



(3.13) 



being Sj — > So for J — > 0. Let us also stress that L'j and L' are non-degenerate, the related 
Hessians being non-vanishing. 

In order to see the main differences between the case here studied and the pure gravity 
one, we have also to consider the energy functions 

m + 1 



E'j 
E' 



-6 exp 



2nm 



v 1 uv + W exp (3f ) + 8irD . 



-6 exp (2v)itv — 3/C exp (i>) + W exp (3v) + 8ttD 



(3.14) 
(3.15) 



We again find that E'j — > E' Q for J — > and /C = 0, once we note that the presence of /C 7^ 
is indeed evident only in the expressions of L and El and their transformed forms. (When 
K 7^ 0, we in fact have to add the term — 3/Cexp (u) in E' .) 

In the following, we will treat separately the two cases with J = and J ^ 0. This 
should make it clearer, on one side, the contribution of a /C-term and, on the other side, the 
range of possibilities tied to choosing different values of fi (and, therefore, of J or n). This 
might also lead to a better understanding of which are the more realistic values that the 
parameter /1 should preferably assume from the physical point of view, also accounting for 
its nature and origin as a Lagrange multiplier. 



A. J = and generic /C 

From Eqs. (13331) and (3.15) (being E' = 0), we find 

v = v{t) = ^hi^t + 2cX (3.16) 

where C\ is an arbitrary integration constant, and 

. , 2 v / 31V(S t + 6d) 5 / 2 - 30^(^0^ + 6d) 3 / 2 + 360nDJ: t + 45C 2 £ 2 . 1 
n = n(t) = ——J , (3.17) 

402^0 

C2 being a further arbitrary integration constant. From now on, in order to show simplified 
formulae, we choose to set C 2 = 0, hence losing the full generality of our solution. (Such a 
generality is indeed guaranteed from the existence of the three integration constants C\, C2, 
and So, plus the arbitrary constant D, which can be seen as a first integral of the condition 
K = 0.) 

Assuming that E is positive, we can rescale time by defining r = S t + 6C1, so that we 
find (with C 2 = 0) 

, v _ 2[r^\Wr - 15/C) + 60y / 37rD(r - 6CQ] 



and 

v(t) = , 
2 V3 



v ( T ) = h n (l). (3.19) 
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From Eqs. (3.7) and (3.8), we thus find 



2t 1 I\U?,t-/ 2 {Wt - 15/C) + 180vr J D(r - 6CO 
a = air) = y - — -= , (3.20) 



2r-V 4 , U/3t 3 / 2 (Wt - 15/C) + 180ttD(t - 6d) 
G — G(t) — y - — = , (3.21) 



so that 



A = A( G (x)) = 1yG -V)= , r , 15 fW /2 . (3.22) 

V V V 7 4[V3WY 5 / 2 - 15^x3/2 + 180tlD(t - 6Ci)] V ' 

On the other hand, the Hubble parameter turns out to be 

H = Hit) = ^1 = H(t) = ^- = 3 ^^ t3/2 ( Wt - 10/C ) + 907rZ ^ ~ 2g i)] f3 23) 
a(t) ' a dr 2t[V3t^ 2 (Wt - 105/C) + 180tiD(t - §C X )\ ' 

In the case of spatially flat universe, Eqs. (3.20), (3.21), (13.221) and (13. 23f) have to be 
rewritten with K = 0, of course. In this case, the expressions so found can also be deduced 
from those we are going to write in the next subsection, once we set J = therein. 

For the time being, let us comment a little on what we have found above, when J = 
and /C 7^ 0. Eq. (3.20) can be more conveniently written as 



a(r) = A^r 3 - Br 2 + Cr^l 2 , (3.24) 

with 

2 fW „ 15/C ^ 180ttL) 
so that the physically meaningful parameters are 

K = -A 2 BY, Q \ W=^A 2 E \ D = ^-A 2 Ci: 2 . (3.26) 

4 4 l07T 



On choosing t = r(t ) = 1 (being t the present time), we get a = a(to) = Ayl — B + C 
and 



,( T ) S *-1 = — -1, ir(r) = £ ^ +2 ^ 3r - 2 ^, (3.27) 

so that the present value of the Hubble term is H = S (6 - 45 + 3C)/(4 - 45 + 4C). 
On the other hand, we can write Eq. (3.21) as 

.4 



G(t) = -^Vr 2 - Bt + Ct 1 / 2 , (3.28) 
v3 



from which Go = G{r = 1) = A a/1 — B + C / \/3. (Note that such a value is in general 
different from the Newtonian one.) With D from Eq. (3.26), we also have the following 
expression for the energy density: 

P= = 16Att(t3 - fir 2 + Ct3/ 2 )3/2 • ( 3 - 29 ) 

Its present value is therefore p — V / 3CS 2 /[16A(1 — B + C) 3 ^ 2 tt). As to the present matter 
content of the universe, we can continue to use the above defined value of Go, and find that 

ft, = ^ = ggii^±a (3.30) 



3# 2 ~3(6-4B + 3C) 2 
only depends on B and C (i.e., not on E ). Of course, it is possible to adjust such parameters 
to get a value of Qq not very far from 0.3. In general, it is in fact possible to see that there 
exist suitable values for B and C giving something not substantially different from what is 
usually produced in the ACDM model, for instance. This will be anyway better outlined in 
the next section, for the case with /C = 0, where a little more refined comment on such a 
possibility is worked out. 

B. J / and K = 



On setting again Ci = and r = £jt + 6Ci, we now have 



v = v(t) 



2n 



(6ra - 1) 



In 



Qn — 1 
12n 



-T 



and 



u = u(t) 



6n 12n — 1 
12— fcl 6n-l ^ g7rjD / _ Q C 
- — 1 t 6n-l -I 



(3.31) 



(3.32) 



(12n-l)E/ Sj 
From them, with C% = we get r = Sjt and a(0) = (which indeed requires attention; see 



Ref. 



291 ] for detailed considerations on that), so that we have 



a = a(r) 



, 6n 12n 

A Bt^- 1 + t 6 ™- 1 



G = G(t) 

where we have defined the constants 



C [r 2 + Bt^, 



3n-l 



.4 
B 

C 



n(l + 6n) 12n 2 



A(n, W, Sj) = 12-^^nT^(6n - 1)~ (12n - l)- n W n ^j- 2n . 



B(n,W,D) = W~ l 2^=^ (3n)^ (Qn - l)^r(12n - l)nD 



C(n, W, Ej) = (6n - l) 2 ^- 1 ) [l2n 2 (12n - 1)] 



1_3n p^ 3n_1 



2(l-3n) 



(3.33) 
(3.34) 

(3.35) 
(3.36) 
(3.37) 
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Thus, we can see that such expressions are generalizations of the ones found in the case 
with J = 0, with the usual caution as to the terms with /C / 0. Also, when D = we 



recover the same results obtained in Ref. 



la ] for the pure gravity model. Of course, let 



us note that above we have had to choose appropriate values of the n parameter, setting 
—1/6, 0, 1/12, 1/6, 1/3. This affects the values of the interaction \x parameter occurring 
in the Lagrangian; we must then assume from the very beginning that our treatment is 
performed with fi ^ 0, 2/27, 2/3. (The values fi = 0, 2/3 have been already ruled out, since 
J ^1,3/2.) 

Furthermore, let us explicitly exhibit the function A = A(r) in the matter-dominated 
gravity regime of the universe when (=^ n 7^ 1/2 and // ^ 6) and /C = 



A = A(r; n) = WG^ = WC^ [t z + Bt^J . (3.38) 

(When J = (=>• n = 1/2) and /C = 0, this reduces to the A-term found in that case.) 
Now the Hubble parameter has the more general expression 

, v ait) u 1 < 1 £ 7 aa Hi + Hot 6 *- 1 

H = Hit) ^AA = n - + -u = H(t) = JL = 1 1 , 3.39 

a(t) u 2 a dr £ jr [#3 + H^—} 



where 



2(9n-l) 

F a = Hi(D) = 48irD(12n - l)n"fe=r , (3.40) 



# 2 = f 2 (^,Ej) Eiy n 



(6n - l) 12n_1 „ 



Yin 



12 

6n 



(3.41) 



# 3 = H 3 (D) = 8ttD[1 + 18n(4n - l)]n^r , (3.42) 
H 4 = H±(W, Ej) = W [S/(6n - 1)] ^ (l2^ - 12^n + S^^^n 2 ) .(3.43) 

In the next section we will use such formulae to produce a cosmological model able to 
be compared successfully with the concordance ACDM one, which makes it interesting to 
further study our model in the future. 



IV. THE K = MODEL AND THE OBSERVATIONS 

A way to understand whether the cosmological model resulting from the above expressions 
for a, G, and A may be considered for further analysis is to study how the results obtained 
can fit observational data from SNIa. Usually one should take at least the Gold Data Set [3^] 

11 



as observational result and choose the constants involved above in order to obtain a best fit. 
Here, we simply limit ourselves to directly compare our K, = model with the concordance 
ACDM one, which, as known, fits those data very well. This should make it possible to infer 
a first feasibility of the model here considered. Of course, further investigations are anyway 
necessary in order to assess all related issues and have to be postponed to future work. 



A. Some useful cosmological parameters 

First of all, let us limit our attention to the /C = case with C-i = 0, so that we can use 
the generic J formulae above. (Things sensibly change when K, ^ and have to be analyzed 
separately.) In order to further simplify the expressions below, let us also set £j = 1 and 
C\ = in them, which arbitrarily fixes the time scale and origin. From r = Ejt = t we then 
have a(0) = 0. Setting W = 1 requires more attention, on the other hand, even if the W 
parameter only appears as a factor in the expression of u. This is, in fact, a very peculiar 
situation which will be dismissed in the last part of the next subsection B, in order to better 
understand how well the model defines the time evolution of G. Here, we can proceed rather 
roughly, because the aim of this subsection is only that of showing that there exists at least 
one reasonable choice of the parameters. Thus, in the following let us start by considering 
the cosmic scale factor a = a(t) as given by Eq. (3.24) with 

n(l + 6n) 12n 2 , , 12n 2 , . , 

A = l2^^n—(Qn-l) — (12n-iy n , (4.1) 

3(l-10n) 6n . . l-12n , , 

B = 2^=^(3n) — (&n-l)-**=r(12n-l)irD, (4.2) 

so that fixing the value of n ^ —1/6, 0, 1/12, 1/6, 1/3 simply yields A as a number, while B 
still depends on D. On the other hand, the expression for G is given by Eq. (3.27) with 

C = (6n - 1) 2(3 "- 1} [I2n 2 (12n - 1)] ^ 3n , (4.3) 

from which, using Eq. (I3.38p . one can get the related expression of A = A(t). 
Other relevant quantities for observations are 

6n 

Hi + Hot 6 "- 1 

H = H(t) = , 1+ 2 6n x , (4.4) 



/ 6n\ 

t iH 3 + H 4 t~j 



O O m = 8nGD - 327rDmi(t)e 1 2^M _ gJV + nj 2 k^ 

ll m — ll m (t) — — — 6n „ 12n • V^.Oj 
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Here, H 1 and H 3 are those in Eqs. (3.33) and (3.35), while H 2 and H4 are now 

"(6ra- l) 12 ™- 1 ' 



# 2 = H 2 (W = Sj = 1) =n 2 



12 



2(3n-l) -1 



#4 = # 4 (W^ = Sj = 1) = (6n - l)^ 1 (^12^t _ I2^r n + 3^^4^n 

On the other hand, the constants fim*- 1 -*, . . . , are all depending on D, while Q r 
not depend on it; they are defined as 



n (2) 



8 2(9n-l) 
-6n 3 1— 6n 

5(l + 6n) 2(12n-l) 1 2n ~ 1 



Q m (1) = 2^3^6^ (1 - 12n) 2 (6n- l)n^^7r 2 L), 



2-l^6^3~ 

8 6re 



2(9n-l) 

(6n — 1) 7rD , 



■ 1 ■. Y-|-l»n inn , n// 2(1 — 9n) 

nJ A > = 2^^3^r(6n- 1)6^ 1 + 2 • 3^^n(4n - 1) 



7TD, 



/-x 4(12n+l) lgn 6n 2(12n-l) 

iZ m l o; = 2 1-6" 3 1 -6«n 1 -6«(6n — 1) e™- 1 



We can also introduce the parameters 

A 



n 



G 



Anu(t)e 2n 



3tf 2 3(2nw + m;) 2 fiyV (3) + n A W^ + Q^a^r 



-2rw + uv) 2 (fi G (1) + ft G (2) ^) 5 



6tf 2 G 2 (2nu + uvf n V) + fi r (4) t^T + fi G (5) t^ 



where 



^A (5) 



[2 18n - 1 (3n) 6n (6n - 1)] 6 "" 1 [1 + 18(4ra - 1)] VD , 



1 2(12n-l) 

-(2n- l)(6n- 1)~bs=t- ? 

2(12ra-l) 

2 (3n) 6 



^ 2 n 2 

7T U 



r(l + 12n) 

[2 3(10n - 1) (3n) 2(9n - 1) (6n - l) 12 "" 1 ]^ (12n - 1)ttD, 



3n(6n — 1) e«- 
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(only and being independent of D) and 

'2 4 3 12n \ 1=6^ 



with, again, only fie*" 2 



n 



6ra 



6 



18n 



[1 + 3n(12n - 5)]tt£) , 
[1 - 18n(l + 6n(2n - 1))] , 



2(6n- 1) 

j- 2 3(7-12n) 3 2(l-3n) n 2(l-12n)] deS ^ _ 12n )V/) 2 , 
l 

2 7-18no2 2(l-9n)l T^6n 

(6n _ 1)6 » J [l + 18n(4n-l)]7rD, 
= [l2 1+3n (6n- l) 1 - 12 "]^^, 

and independent of D. 



(4.20) 

(4.21) 
(4.22) 
(4.23) 
(4.24) 



B. Comparison with the ACDM model 

In order to begin to compare the model characterized by the parameters above with 
the usual ACDM model, let us consider the behaviour of the scale factor a for large and 
small times. In the former case, this is equivalent to take the expression in Eq. (3.26) with 
B — > 0, a t ^oo = At6^, while in the latter we can assume a t ^ = AB n t^^ . This means 
that the two exponents characterizing the asymptotic time behaviours of the scale factor are 
respectively 



Yin 1 



6n z 



where we want that p\ > 1 and p 2 < 1, which implies that we find a limited range of 
variability for the n parameter, (3 — v / 3)/6 < n < (3 + \/3)/6. This gives rise to a constraint 
on the values of the fx parameter introduced in the Lagrangian at the beginning of our 
investigation, since fx is indeed equal to 2 at both extremals, but we also find that fx — > oo 
when n — > 1/3. Thus, the range of values of fx is indeed open from above, i.e. fx > 2. 

Let us now consider the behaviour of G — G(t) for large and small times. As before, 
this yields G^ = Ct 2 ^^ (for B -> 0) and G^o = (75 3 "-4^. The two exponents 
involved are now 



p 3 = 2(3n - 1) , p 4 



2(1 -3nf 



(4.26) 



6n — 1 

Here, if we want that both exponents do not attain too high a value, we can stick around 
n — 1/3, which, as already above, represents a degenerate situation and cannot be considered 
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exactly in our analysis. Anyway, let us underline that, for n — > 1/3, we get 
a(t) — A[t 2 (B + t 2 )]^ , -> 2 f , + 4t ,l , G(t)^C. 



(4.27) 



/J ' w 3Et + 3t 3 ' 
which recovers (once we set B = 1) the solution found in Ref. [29j, already known to fit 



observational Snla data 
can be written as 



29 



321 ] . In general, however, H = H(t) is given by Eq. (4.4) and 



H(t) 



6n B + 2t~ 



(4.28) 



(6n - l)t ( -B + t ^ 

Once we choose the present time to = 1, the present scale factor and Hubble term respectively 
become 

a = A(l + B)\ #o = 6( ^t,f )W ' , (4-29) 



% 1 + S)(6n-1) 

so that in such units H ~ 1. In what follows, we therefore choose to set H = 1 exactly, even 
if this has to be considered as simply heuristic, since we have already fixed the time scale 
by setting Ej = 1 and the age of the universe with to = 1. This makes things analytically 
simpler but, at the same time, underlines the fact that we have chosen to consider the model 
only at a first glance, without even trying to fit observational data, just to look at a first 
viability of our model. A more correct analysis then deserves more and better work, of 
course. Here, we have chosen to focus our considerations on the analytical method adopted 
to get the model itself rather than on its compatibility with observations. 
If Hq = 1 one finds that 

B 



■1 + 6n - Yin 2 



1 — 6n + Qn 2 

which inserted into the expression of the Hubble term leads to 

-l+6n-12n 2 



(4.30) 



6rr 



H 



1— 6n+6n 2 



6ra 

+ 2*6^-1 



(6n - l)t 



-l+6n-12n2 + 



(4.31) 



1— 6n+6n 2 

Let us then note that we have not chosen to normalize a , and that, using such a value for 
B, the definition of the redshift z = ao/a — 1 leads to 



6n 2 



6n + 6n 2 



9\ 

-1 + 6n — 12n z )t6«-i 

i 1- £6n-l 

1 — 6n + 6n 2 



(4.32) 



In this way, we see that the only parameter to be considered is n. 
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0.25 0.5 0.75 1 1.25 1.5 



FIG. 1: Comparison between the Hubble terms H = H(z) in our model (red) and H = H{z) in 
the ACDM model (black). 




0.25 0.5 0.75 1 1.25 1.5 



FIG. 2: Comparison betwwen the luminosity distances di = cIl(z) in our model (red) and cIl = 
d L (z) in the ACDM model (black). 

Let us now compare the function H = H(z) we can deduce from such expressions (where 
we set, for example, n = 0.3) with that usually estimated in a ACDM model with Cl m = 0.27 

H A = y / l-fi m + fi m (l + z 2 ) 3 , (4.33) 

hence finding a first but apparently not very good agreement, as shown in Fig. 1. (Notice 
that the tilde parameters are those peculiar to the ACDM model, which can of course be 
different from the ones in our own model.) As a matter of fact, however, using the Hubble- 
free luminosity distance gives a plot (see Fig. 2) with a very good agreement, at least for 
z < 2. 

To understand better what goes on, let us first of all stress that Ej = 1 and t — 1 tell 
us that we are using the age of the universe as unit. On the other hand, if we now dismiss 

16 



our previous choice for W, write the expression of B without setting W = 1, and compare it 
with the one we have found as a result of Hq = 1, we get 1.077 = 2A09.08D/W , from which 
D = 0.00045W. Making something similar with the expression of G analogously leads to 
C = 1.16/W 01 . In this way, it turns out that the present value of G — G(t) is 

Go = G(f = 1) = ^ ■ (4-34) 

Of course, we do not expect that Go = Gn = ^Newton, even if we can determine W in order 
to get it. We can take W = 1 or, for example, G = G^ = 1; in this latter case, we find 
W = 2.11. Anyway, we have to consider that small differences between Go and G^ make 
this value vary a lot. To better investigate this issue, we can evaluate the fractional time 
rate of change of G 

G 

— = -8.33 x 1(T 2 . (4.35) 
G 

Here, it is important to remember that we are using the age of the universe as unit, so that 
the effective rate is of the order 10~ 13 yr -1 . We can also estimate how G (with n = 0.3) 
varies from now (z = 0) up to the equivalence time (z eq ~ 10000), by exploiting the formula 



G(t) = TWT-7-2 ■ ( 4 -36) 

jO.25 ~T L 

This shows that the variation is 30% and can be seen as acceptable in such first rough 
considerations. On the other hand, the infrared fixed-point hypothesis is that the nontrivial 
running is due to quantum fluctuations with momenta smaller than fc cosmo , where l/fc COS mo 
is the length scale characterizing the largest localized structures in the Universe. There is 



no conflict with the classical tests of general relativity [31]. 

Looking at the present value of A, we find Ao = A (to = 1) = 2.11, which is indeed in Hq 2 
units. As to the more fundamental test on Q m , we have 



__ 8tiG D _ 2.09 x 10- 3 H/°- 9 
fim0 = ~ A* • (437) 

Since A = 2.07 x 10" 1 x W - 93 , we get 

tt m0 = 2.33 x 10" 1 , (4.38) 
which is indeed a very reasonable value, because all our comparison procedure was meant 



to emulate what done with the solution obtained in Refs. 
does not depend on W .) 

17 



29 



32]. (Notice that this result 



We can also evaluate 

fl A0 = = 7.04 x 10" 1 , Q G0 = ^—.{t = t ) = 6.25 x 10~ 2 , (4.39) 
3 oG z 

from which we deduce the usually expected constraint 

O m0 + ^ao + fico = 1 • (4.40) 



V. CONCLUSIONS 



Since a renormalization-group approach seems to show that quantum Einstein gravity 



is asymptotically safe 
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despite being non-renormalizable at perturbative level, we began 
investigating the cosmological applications in Ref. [ijj], devoted to the Lagrangian and 
Hamiltonian formulation with variable G and A. In the present paper, we have solved 
the equations from an RG-improved gravity Lagrangian for a homogeneous and isotropic, 
matter-dominated universe, finding the related dynamical behaviour for the gravitational 
coupling G and the cosmological A-term. We have used the so-called Noether Symmetry 
Approach, yielding a coordinate transformation which leads to a form of the equations easily 



and exactly solvable [l8|, Q, This has required to impose the existence of a Noether 
symmetry for the point-like Lagrangian describing the cosmological dynamics. The existence 
of the coordinate transformation so found, which might indeed also be guessed a priori, is 
a direct product of the machinery worked out here even if the solutions we obtain are 
independent of the very existence of the symmetry. On the other hand, as we have already 
found in the pure gravity case [18], such a symmetry brings into our treatment some strong 
conditions, like those on the values of K, and the J parameter, and on the form of the 
function A = A(G). 

While in Ref. [17] a power-law behaviour for the scale factor was guessed from the 
beginning to solve the equations, we have seen already in Ref. [l8j] that this can instead 
result exactly and generally from the method itself. When K, = 0, in Ref. [l7| it was assumed 



a = At a for the scale factor; with arbitrary A, this gives a± = (3 ± y/9 + 12£ 2 A*)/6, in close 
connection with the hypothesis of being in the neighbourhood of a fixed point (g+, A*), which 
constrains G and A to be 

G = G(t) = g^"H 2 , A = kit) = \^H' 2 GA = #*A* = constant . (5.1) 
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This allows arbitrarily large values for a+ (since £ is undetermined), and constrains the 
interaction parameter to be /x± = 3a±/2 = (3 ± a/9 + 12£ 2 A*)/4 17(. If A* > 0, this leads 
to power-law inflation for the "+" solution. 



Under such assumptions it is possible, as in Ref. [18|, to see that the general asymptotic 
trend of a for large r is power-law, a ~ r p (with p = (3 — 2J) 2 /[3(J 2 — 3 J + 2)]), without 
strictly imposing to be near a fixed point, but rather being probably very far from it. We 
again find acceleration only when J<lorJ>2,p being always > 1 in these ranges; when 
J is near the values 1 and 2, the exponent can assume any large value, therefore yielding a 
possible strong power-law behaviour. On the other hand, the interaction parameter /x is a 
function of J such that both p and /i assume symmetric values for J < 1 and J > 2. We 
also find a power-law behaviour, with exponent p' = 1/(1 — J), for the function G = G(r). 
A special case is obtained when p' = 2, that is J = 1/2; here, we have an accelerated stage, 
with p = 16/9. In the J = case, p is instead fixed to be 3/2, with a soft acceleration for 
the universe. 

It is interesting to note that the model with K, ^ also makes it possible to perform 
a sort of rough comparison with the ACDM model, even if more refined work has to be 
done in such a direction. Here, in fact, we have limited ourselves to outline the procedure 
adopted to find solutions rather than the necessary comparison with observation, which is 
then postponed to future work. 

In general, we have to stress that the presence of an ordinary matter (dust) term in the 
theory makes larger differences with respect to what is found in the proximity of the non- 
Gaussian ultraviolet fixed point, and also with respect to the pure gravity case (which can 
be probably interpreted as characterizing the epoch soon outside the attraction basin of that 
fixed point), and can be a little better understood by means of the solutions resulting from 
the cosmological equations. As a matter of fact, we get power- law dependence on time for the 
scale factor, which might also allow for accelerated trends. Of course, these would be now set 
up in a period next to the much earlier one characterized by a pure gravity regime, in which 
standard inflation was born [18] . This could mean that what we have effectively treated here 
regards a cosmic evolutionary stage close to the current one, this latter being first assumed 
(and tested) as matter-dominated and then recently discovered as an accelerated one. In 
our context, quantum effects induce running of both the gravitational coupling G = G(t) 
and the cosmological term A = A(t). This in fact yields primordial inflation soon after the 
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universe exits the region where the attraction basin of the non-Gaussian (ultraviolet) fixed 
point works [18j, as well as another inflationary behaviour in or, probably, after a matter- 
dominated epoch, and always without need to introduce a scalar field in the cosmic content. 
The usual matter domination should indeed occur in between, also after what is usually 
described as radiation era. This means that we do not have any concrete continuous model 
of the evolution of the universe, while also assuming that our considerations might be valid 
once we accept an onset of quantum effects only at the beginning and at the end of such 
evolution. 

What also lacks in this sort of a patchwork reconstruction of cosmological history mainly 
concerns the way one could establish suitable links among these separated stages. Nonethe- 
less, if we do accept the very recently assumed existence of another non-Gaussian, infrared 
fixed point late in the history of the universe, it could indeed characterize a later inflationary 
period. (This last regime of the evolution of the universe, typically pictured as an asymp- 
totic dark-energy dominated stage, is mostly described by means of a scalar field. Let us 
again stress, here, that in our treatment we do not need to introduce such an ingredient 
in the cosmic content.) Thus, other important physics would be needed, linking this later 
era to the previous ones, which prevents a thorough understanding of the whole picture we 
are drawing. As a result, however, we would be partially facing the construction of a global 
patchwork model for the cosmic evolution which seems to be of some interest, because, once 
we introduce the quantum effects in a RG-inspired effective theory of gravity, we do not 
need any other mechanism to work out accelerated stages. 

Still lacking is also the concrete possibility to generate an acceptable radiation-dominated 
regime, which is instead needed to describe what we understand as the usual standard 
cosmological model, in the framework of RG-inspired cosmology. The procedure adopted 
in this paper, using the Noether Symmetry Approach, does not in fact easily work with a 
Lagrangian where the matter term is L m = Da^ 1 (since 7 = 4/3 for radiation). Thus, other 
methods have to be worked out for treating the problem of solving cosmological equations 
in such a period. Finally, it should be stressed that, in the context of our paper, the 
development of a perturbation theory able to produce suitable structures in the universe is 
meaningless, since it can indeed work only in an appropriate matter-dominated era, when 
the scale factor is a ~ t 2 / 3 , which we never find, instead, in our treatment. 
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